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Abstract. Spin-glass and chiral-glass orderings in three- dimensionaH 
Heisenberg spin glasses are studied both by equilibrium and off- 
equilibrium Monte Carlo simulations. Fully isotropic model is 
found to exhibit a finite-temperature chiral-glass transition with- 
out the conventional spin-glass order. Although chirality is an 
Ising-like quantity from symmetry, universality class of the chiral- 
glass transition appears to be different from that of the standard 
Ising spin glass. In the off-equilibrium simulation, while the spin 
autocorrelation exhibits only an interrupted aging, the chirality 
autocorrelation persists to exhibit a pronounced aging effect rem- 
iniscent of the one observed in the mean-field model. Effects of 
random magnetic anisotropy is also studied by the off-equilibrium 
simulation, in which asymptotic mixing of the spin and the chi- 
rality is observed. 



§1. Introduction 

Ordering of complex systems has attracted interest of re- 
searchers working in the field of numerical simulations. Well- 
known examples of such complex systems may be a variety of 
'glassy' systems including window glasses, orientational glasses of 
molecular crystals, vortex glasses in superconductors and spin- 
glass magnets. Often, in the dynamics of such complex systems, 
characteristic slow relaxation is known to occur. It has been a 
great challenge for researchers to clarify the nature and the origin 
of these slow dynamics, as well as to get fully equilibrium proper- 
ties by overcoming the slow relaxation. In particular, spin glasses 
are the most extensively studied typical model system, for which 
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numerous analytical, numerical and experimental works have been 
made [1]. 

Studies on spin glasses now have more than twenty years of 
history. Main focus of earlier studies was put on obtaining the 
equilibrium properties of the spin-glass ordering. Owing to ex- 
tensive experimental studies, it now seems well-established that 
the spin-glass magnets exhibit an equilibrium phase transition at 
a finite temperature [1]. From theoretical side, there now seems 
to be a consensus that the lower critical dimension (LCD) of an 
Ising spin glass is between d — 2 and 3, while the LCD of vector 
spin glasses is greater than d = 3. In other words, at d = 3, only 
an anisotropic Ising spin glass exhibits an equilibrium spin-glass 
transition at a finite temperature [2,3], whereas isotropic Heisen- 
berg spin glass exhibits only a zero-temperature transition [4-7]. 
Meanwhile, it has been known that the magnetic interactions in 
many of real spin-glass materials are Heisenberg-like, in the sense 
that the magnetic anisotropy is much weaker than the exchange 
energy. Apparently, there is a puzzle here: How can one reconcile 
the absence of the spin-glass order in an isotropic Heisenberg spin 
glass with the experimental observation? While weak magnetic 
anisotropy inherent to real materials is often invoked to explain 
this apparent discrepancy, it remains puzzling that no detectable 
sign of Heisenberg-to-Ising crossover has been observed in experi- 
ments which is usually expected to occur if the observed spin-glass 
transition is caused by the weak magnetic anisotropy [1,4]. 

In order to solve this apparent puzzle, a chirality mecha- 
nism of experimentally observed spin-glass transitions was re- 
cently proposed by one of the authors [6], on the assumption 
that an isotropic 3D Heisenberg spin glass exhibited a finite- 
temperature chiral-glass transition without the conventional spin- 
glass order, in which only spin-reflection symmetry was broken 
with preserving spin-rotation symmetry. 'Chirality' is an Ising- 
like multispin variable representing the sense or the handedness 
of the noncoplanar spin structures (more detailed definition will 
be given below). It was argued that, in real spin-glass magnets, 
the spin and the chirality were "mixed" due to the weak magnetic 
anisotropy and the chiral-glass transition was then "revealed" via 
anomaly in experimentally accessible quantities. Meanwhile, theo- 
retical question whether there really occurs such finite-temperaturej 
chiral-glass transition in an isotropic 3D Heisenberg spin glass, a 
crucial assumption of the chirality mechanism, remains somewhat 
inconclusive [6,7]. 

More recently, there arose a growing interest both theoreti- 
cally and experimentally in the off- equilibrium dynamical prop- 
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erties of spin glasses. In particular, aging phenomena observed 
in many spin glasses [8] have attracted attention of researchers 
[9,10]. Unlike systems in thermal equilibrium, relaxation of phys- 
ical quantities depends not only on the observation time t but 
also on the waiting time tyj, i-e., how long one waits at a given 
state before the measurements. Recent studies have revealed 
that the off-equilibrium dynamics in the spin-glass state generally 
has two characteristic time regimes [9-11]. One is a short-time 
regime, to << t « (to is a microscopic time scale), called 
'quasi-equilibrium regime', and the other is a long-time regime, 
t » t,^, called 'aging regime' or 'out-of-equilibrium regime'. In 
the quasi-cquilibrium regime, the relaxation is stationary and the 
fluctuation-dissipation theorem (FDT) holds. The autocorrelation 
function at times and t + tyj is expected to behave as 

C(^„,^ + U^g^^ + ^^g^^ (i) 

where q^^ is the equilibrium Edwards- Anderson order parameter. 
In the aging regime, the relaxation becomes non-stationary, FDT 
broken, and the autocorrelation function decays to zero as t ^ oo 
for fixed t^. 

On theoretical side, both analytical and numerical studies 
of off-equilibrium dynamics of spin glasses have so far been lim- 
ited to Ising-like models, including the Edwards-Anderson (EA) 
model with short-range interaction [12-14] or the mean-field mod- 
els with long-range interaction [11,15-17]. Although these anal- 
yses on Ising-like models succeeded in reproducing some of the 
features of experimental results, it is clearly desirable to study 
the dynamical properties of Heisenberg-like spin-glass models to 
make a direct link between theory and experiment. 

In the present article, we report on our recent results of 
equilibrium as well as off-equilibrium Monte Carlo simulations 
on isotropic and anisotropic 3D Heisenberg spin glasses. Order- 
ing properties of both the spin and the chirality will be studied, 
aimed at testing the validity of the proposed chirality scenario 
of spin-glass transitions. We note that Monte Carlo simulation is 
particularly suited to this purpose, since, at the moment, chirality 
itself is not directly measurable experimentally. By contrast, in 
numerical simulations, it is quite straightforward to measure the 
chirality. 

§2. Chirality 
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Frustration in vector spin systems such as the XY and Heisen-| 

berg models often causes noncollinear or noncoplanar spin struc- 
tures. Such noncolhnear or noncoplanar orderings give rise to a 
nontrivial chirality. In the case of two-component XY spins, one 
may define a local chirality for the two neighboring spins at the 
i- and j-th sites by Kij = x Sj \z— S^Sj — S^Sj. When the 
spin configuration in the ordered state is noncollinear, the local 
chirality Kij defined above takes a nonzero value. 

In the case of three-component Heisenberg spins, three spins 
are necessary to define a scalar chirality. Thus, in the Heisenberg 
case, one may define a local scalar chirality for three neighboring 
spins (spin triad) at the i-, j- and k-th sites by Xijk — Si • Sj x Sfe. 
It takes a nonzero value for any noncoplanar spin configurations 
but vanishes for any planar spin configurations. Note that, in 
either case, the chirality is a pseudoscalar in the sense that it is 
invariant under global spin rotation but changes sign under global 
spin refiection. Possible chiral ordering is related with a breaking 
of the refiection symmetry with preserving the rotation symmetry. 

The model we simulate is the classical Heisenberg model on 
a simple cubic lattice with the nearest-neighbor random Gaussian 
couplings, Jij and D^^^ , defined by the Hamiltonian 

7i = - J](J,,Si-S,+I?^^5f5r), (2) 

<ij> 

where ={Sf , Sf , Sf) is a three-component unit vector, and the 
sum runs over all nearest-neighbor pairs with N — LxLxL spins. 
Jij is the isotropic exchange coupling with zero mean and variance 
J, while -D^*^ {fJ,^ = x,y,z) is the random magnetic anisotropy 

with zero mean and variance D which is assumed to be symmetric 
and traceless, L>f/ = L*,^ and L'f/ = 0. 

Wc define the local chirality at the i-th site and in the fi-th 
direction, x^^, for three Heisenberg spins by, 

Xin — Si+e^ • (Sj X Si_e^), (3) 

where {/i = x, y, z) denotes a unit lattice vector along the fx- 
axis. 



§3. Equilibrium simulations 
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First, we report on our equilibrium Monte Carlo simulation of 
a fully isotropic 3D Heisenberg spin glass defined by eq.(2) with 
D = 0. Monte Carlo simulation is performed based on an 'ex- 
tended ensemble' method recently developed by Hukushima and 
Nemoto [18], where the whole configurations at two neighboring 
temperatures of the same sample are occasionally exchanged with 
the system remaining at equilibrium. By this method, we suc- 
ceeded in equilibrating the system down to the temperature con- 
siderably lower than those attained in the previous simulations. 
We run in parallel two independent replicas with the same bond 
realization and compute an overlap between the chiral variables 
in the two replicas. 



In terms of this chiral overlap, q^, the chiral-glass order parameter. 



q^Q, and the Binder cumulant of the chirality, gcG, SiXe calculated 

by 



where < ■ ■ ■ > represents the thermal average and [■ ■ ■] represents 
the average over bond disorder. At the possible chiral-glass tran- 
sition point, curves of gcG against T for different L should merge 
or cross asymptotically for large L. 

For the Heisenberg spin, one can introduce an appropriate 
Binder cumulant in terms of a tensor overlap q^^, {fi, v = x,y, z) 
which has 3^ = 9 independent components. 





= [< 4 >]^ 



(5) 






(7) 



via the relation. 



^SG = t; 




(8) 
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The lattice sizes studied are L = 6, 8, 10, 12, 16 with periodic 
boundary conditions. In the case of L = 12, for example, we 
prepare 50 temperature points distributed in the range [0.08 J, 
0.25J] for a given sample, and perform 4.7 x 10^ exchanges per 
temperature of the whole lattices combined with the same number 
of standard single-spin-flip heat-bath sweeps. For L = 12, we 
equilibrate the system down to the temperature T/J = 0.08 J, 
which is lower than the minimum temperature attained previously. 
Sample average is taken over 1500 (L = 6), 1200 (L = 8), 640 
(L = 10), 296 (L = 12) and 32 (L = 16) independent bond 
realizations. Equilibration is checked by monitoring the stability 
of the results against at least three-times longer runs for a subset 
of samples. 

The size and temperature dependence of the Binder cumu- 
lants of the spin and of the chirality, ^fsG and gcG, are shown in 
Fig. 1(a) and (b), respectively. As can be seen from Fig. 1(a), ^rsG 
constantly decreases with increasing L at all temperatures stud- 
ied, suggesting that the conventional spin-glass order occurs only 
at zero temperature, consistent with the previous results [1,4-7]. 
A closer inspection of Fig. 1(a), however, reveals that qsg for larger 
lattices (L = 10, 12, 16) exhibits an anomalous "upturn" around 
T/J ~ 0.1 — 0.15, suggesting that a change in the ordering behav- 
ior occurs in this temperature range. 
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FIG.l Temperature and size dependence of the Binder cumulants of 
the spin (a) and of the chirality (b) of a 3D isotropic Heisen- 
berg spin glass. 

As can be seen from Fig. 1(b), the curves of gcc for differ- 
ent L do not cross, but show a tendency to merge for larger L in 
the temperature range where the curves of ^rsG exhibit an anoma- 
lous upturn. However, it is not possible to determine only from 
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the data of ^fcG whether the system exhibits a finite-temperature 
transition into the chiral ordered state with some "critical" charac- 
ter, or it exhibits only a zero-temperature transition with rapidly 
growing correlation length. 

Thus, we also have tried the standard finite-size scaling anal- 

ysis for the chiral-glass order parameter Qqq, by adjusting the 
transition temperature Tqg and the exponents PcG and ucG as 
fitting parameters. The best fit is obtained for Tcg/J ~ 0.176, 
vcG ~ 1-15 and /Jcg/'^cg ~ 1.30, with the associated x^^f-value, 
Xdof ~ ^-3; see Fig. 2(a). If, on the other hand, the same data are 
fitted to the finite-size scaling form expected for a T = transition 
with nondegenerate ground state, i.e.. Tec = and /^cc/i^CG = 0, 
the best fit is obtained for vcg ~ 3.2, but with the associated 
X^-value, ~ 29.5, which is significantly larger than the best 
X^-value obtained with assuming Tqg > 0. Thus, while our indi- 
vidual estimates of Tcg, '^cg and /Jcg/'^cg might not be so accu- 
rate since we are trying here to determine the three intercorrelated 
fitting parameters simultaneously, our analysis strongly favors a 
finite-temperature chiral-glass transition over a zero-temperature 
transition. Indeed, this conclusion is corroborated with the results 
of the dynamical simulation to be presented in the next section. 




FIG. 2 Finite-size scaling plots of the chiral-glass order parameter of 
a 3D isotropic Heisenberg spin glass, with assuming a finite- 
temperature transition Tqg > 0. The chiral-glass order pa- 
rameter is divided by the magnitude of the local chirality.[6,7] 

In Fig. 3, we display the distribution function of the chiral- 
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overlap defined by 



Pcoiqx) = [< ^iQx - ^x) >i (9) 

calculated at a temperature T/J = 0.1, below the estimated 
chiral-glass transition temperature. The shape of the calculated 
Pcg{Qx) somewhat different from the one observed in the stan- 
dard Ising-like models such as the 3D EA model or the mean-field 
SK mode. Pcg(<?x) standard 'side-peaks' corresponding to 
the Edwards- Anderson order parameter ig^Q, which grow and 
sharpen with increasing L as is usually the case in the spin-glass 
state. In addition to the side peaks, a 'central peak' at = 
shows up for larger L, which also grows and sharpens with in- 
creasing L. This latter aspect, i.e., the existence of a central peak 
growing and sharpening with the system size, is a peculiar feature 
of the chiral-glass ordered phase never observed in the EA model 
or the SK model: It might suggest a nontrivial structure in the 
phase space associated with the chirality. This peculiar feature 
is reminiscent of the behavior characteristic of some mean-field 
models showing the so-called one-step replica- symmetry breaking 
(RSB). Thus, it is tempting to deduce that the chiral-glass phase 
of a 3D Heisenberg spin glass has a character of such one-step 
RSB. If the ordering is really of such type, no crossing of the 
Binder ratio needs to occur at T = Tea, and our data of ^CG are 
consistent with the scenario. Further studies are in progress to 
clarify the situation. 

0.05 I < < < < < < < < — I 




FIG. 3 Chiral-overlap distribution function of a 3D isotropic Heisen- 
berg spin glass below Tec- The temperature is T/J = 0.1. 
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§4. OflF-equilibrium simulations 

In this section, we report on the results of off- equilibrium 
Monte Carlo simulations both on isotropic {D = 0) and anisotropic! 
[D > 0) models. Unlike the case of equilibrium simulation, the 
system here is never in full thermal equilibrium. Recent studies 
have revealed that one can still get many useful information from 
such off-equilibrium simulations, even including certain equilib- 
rium properties. The quantities we are mainly interested here are 
the spin and chirality autocorrelation functions defined by 

Csitn;, t + t^) = ^ ^[< Si{t^) ■ Si{t + t^) >], (10) 



C^{tw,t + ty,) = ^X![< X^^^{tw)X^^^{t + t^,) >]. (11) 

Monte Carlo simulation is performed based on the standard 
single spin-flip heat-bath method. Starting from completely ran- 
dom initial configurations, the system is quenched to a working 
temperature. Total of about 3 x 10^ Monte Carlo steps per spin 
[MCS] are generated in each run. Sample average is taken over 30- 
120 independent bond realizations, four independent runs being 
made using different spin initial conditions and different sequences 
of random numbers for each sample. The lattice size mainly stud- 
ied is L = 16 with periodic boundary conditions, while in some 
cases lattices with L — 12 and 24 are also studied. 

Let us begin with the fully isotropic case, D = 0. The 
spin and chirality autocorrelation functions at a low temperature 
T/J = 0.05 are shown in Fig. 4 as a function of t. For larger t^,, 
the curves of the spin autocorrelation function Cs come on top of 
each other in the long-time regime, indicating that the stationary 
relaxation is recovered and aging is interrupted. This behavior 
has been expected because the 3D Heisenberg spin glass has no 
standard spin-glass order [1,4-7]. Similar interrupted aging was 
observed in the 2D Ising spin glass which did not have an equilib- 
rium spin-glass order [13]. By contrast, the chiral autocorrelation 
function shows an entirely different behavior: Following the 
initial decay, it exhibits a clear plateau around t ~ t^, and then 
drops sharply for t > It also shows an eminent aging effect, 
namely, as one waits longer, the relaxation becomes slower and 
the plateau-like behavior at t ~ t^, becomes more pronounced. 
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FIG. 4 Spin (a) and chirality (b) autocorrelation functions of a 3D 
isotropic Heisenberg spin glass at a temperature T/J = 0.05 
plotted versus log^^g ^ various waiting times tw The lattice 
size is L = 16 averaged over 66 samples. 




FIG. 5 The same data as in Fig. 4, but plotted versus logiQ{t/tyj). 

In Fig. 5, Cs and C-^ are replotted as a function of the scaled 
time t/tw Reflecting its interrupted aging, the curves of Cs for 
larger now lie below the ones for smaller (subaging). By 
contrast, the curves of C-^^ for various cross around t/tyj ~ 1, 
and at t > i^, the data for larger lie above the ones for smaller 
t-uj (superaging). Such super aging behavior of means that the 
aging in chirality is more enhanced than the one expected from the 
naive t/tyj-scaling. Note that, although the chirality is an Ising- 
like variable from symmetry, the observed superaging behavior is 
in contrast to the aging behavior of the 3D EA Ising model which 
was found to satisfy a good t/tw-scalmg in the aging regime [13]. 
It should also be noticed that the plateau-like behavior observed 
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here has been hardly noticeable in simulations of the 3D EA Ising 
model. Rather, the behavior of observed here is reminiscent 
of the one observed in the mean-field model such as the SK model 
[15-17]. This correspondence might suggest that an effective in- 
teraction between the chiralities is long-ranged. Indeed, at least 
in case of XF spins, the interaction between chiralities is known 
to be Coulombic [19], while such analytical information is not 
available for Heisenberg spins. 

While the plateau-like behavior observed in is already 
suggestive of a nonzero chiral Edwards-Anderson order parameter, 
?CG ^ 0' more quantitative analysis similar to the one recently 
done by Parisi et al for the 4D Ising spin glass [14] is performed to 
extract q^Q from the data of in the quasi-equilibrium regime. 

Finiteness of g^Q is also visible in a log-log plot of versus 
t as shown in the inset of Fig. 6, where the data show a clear 
upward curvature. We extract q^Q by fitting the data of for 
tyj — ^ X 10^ to the power-law form of eq.(l) in the time range 
40 < t < 3, 000 satisfying t/t^ < 0.01. Stabihty of the result has 
been checked by examining the robustness of the result against 
the change in the value of t^j and the time window used in the fit. 
The obtained g^G' plotted as a function of temperature in Fig. 6, 
clearly indicates the occurrence of a finite-temperature chiral-glass 
transition at Tea/ J — 0.157 ± 0.01 with the associated order- 
parameter exponent PcG = 1.1 ± 0.1. The size dependence turns 
out to be rather small, although the mean values of Qqq tend 
to slightly increase around Tec with increasing L. Since both 
finite-size effect and finite-^^t, effect tend to underestimate Qqq, 
one may regard the present result as a rather strong evidence of 
the occurrence of a finite-temperature chiral-glass transition. 
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FIG. 6 Temperature dependence of the Edwards-Anderson order pa- 
ll 




ramctcr of the chirality of a 3D isotropic Heisenberg spin 
glass. The data are averaged over 30-120 samples. Inset ex- 
hibits the log-log plot of the t-dependence of the chirality 
autocorrelation function in the quasi-equilibrium regime for 
L = 16 and L„ = 3 x 10^. 



If one compares the present estimate with the static result in 
the preceding section, the estimated transition temperature is a 
bit lower than, but is roughly consistent with the static estimate 
Tcg/J 0.17, while the exponent PcG is somewhat smaller than 
the static estimate, Pqg ~ In any case, the obtained exponent 
PcG ~ 1.1 (or still larger value from the static estimate) is consid- 
erably larger than the value of the 3D EA model P ~ 0.5 [1-3] , and 
is rather close to the value of the mean- field model /? = 1. This 
suggests that the universality class of the chiral-glass transition 
of the 3D Heisenberg spin glass might be different from that of 
the standard 3D Ising spin glass. According to the chirality mech- 
anism, the criticality of real spin-glass transitions should be the 
same as that of the chiral-glass transition of an isotropic Heisen- 
berg spin glass, so long as the magnitude of random anisotropy is 
not too strong. If one tentatively accepts this scenario, the present 
result opens up a new interesting possibility that the universality 
class of many of real spin-glass transitions might differ from that 
of the standard Ising spin glass, contrary to common belief. 

In the presence of weak anisotropy D > 0, chirality scenario 
predicts at the static level that the transition behavior of chirality 
remains essentially the same as in the isotropic case, whereas the 
spin is mixed into the chirality, asymptotically showing the same 
transition behavior as the chirality [6]. In order to see whether 
such "spin-chirality mixing" occurs in the off-equilibrium dynam- 
ics, further dynamical simulations are performed for the models 
with random anisotropies D/J = 0.01 ~ 1. While chirality ex- 
hibits essentially the same dynamical behavior as in the isotropic 
case (not shown here), the behavior of spin ai t > changed 
significantly in the presence of anisotropy. As an example, the 
spin autocorrelation in the case of weak anisotropy D/J= 0.01 
is shown in Fig. 7. Even for such small anisotropy, spin is found 
to show superaging behavior asymptotically at t >> similar to 
that of the chirality in zero and weak anisotropies, demonstrating 
the spin-chirality mixing. 
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FIG. 7 Spin autocorrelation function of the weakly anisotropic 3D 
Heisenberg spin glass with D/ J = 0.01 plotted versus \ogiQ{t/t 
The lattice size is L = 16 averaged over 60 samples and the 
temperature is T/J = 0.05. 
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FIG. 8 Zero- field-cooled magnetization of an anisotropic 3D Heisen- 
berg spin glass with D/J = 0.05 plotted versus log^^g ^- The 
field is h/J — 0.05 and the temperature is T/J — 0.05. The 
lattice size is L = 16 averaged over 80 samples. 

Experimentally, thermoremanent magnetization (TRM) or 
zero-field-cooled (ZFC) magnetization is found to show an ap- 
proximate t/tu,-scaling in the aging regime, with small deviation 
from the perfect scaling in the direction of subaging [9] . Although 
this seems in apparent contrast to the present result, it should 
be noticed that standard aging experiments have been made by 
measuring the magnetic response, not the autocorrelation. Re- 
cent numerical simulation by Yoshino et al revealed that, at least 
in the case of the SK model, TRM showed the subaging even 
when the spin correlation showed the superaging [20]. Thus, wc 
also calculate the ZFC magnetization for an anisotropic model 
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with D / J = 0.05: After the initial quench, the system is evolved 
in zero field during t^j MCS. Then, an external field of intensity 
H/ J = 0.05 is turned on and the subsequent growth of the magne- 
tization M{t] tw) is recorded. As can be seen from Fig. 8, the data 
show the near t/t^-scaling in the aging regime t > tyj where the 
spin-autocorrelation shows the superaging. Thus, the observed 
tendency is roughly consistent with experiments. It might be in- 
teresting to experimentally investigate the aging properties of spin 
correlations of Heisenb erg-like magnets in search for possible su- 
peraging behavior. 

§5. Summary 

In summary, spin-glass and chiral-glass orderings in 3D Heisen-§ 
berg spin glasses are studied with and without random anisotropy 
by Monte Carlo simulations. The results are basically consistent 
with the chirality mechanism: In the isotropic case, clear evi- 
dence of the occurrence of a finite-temperature chiral-glass tran- 
sition without the conventional spin-glass order is presented both 
by equilibrium and off-equilibrium simulations. Spin and chirality 
show very different dynamical behaviors consistent with the 'spin- 
chirality separation'. While the spin autocorrelation exhibits only 
an interrupted aging, the chirality autocorrelation persists to ex- 
hibit a pronounced aging effect reminiscent of the one observed 
in the mean-field model. The universality class of the chiral-glass 
transition appears to be different from that of the the standard 
Ising spin glass. In the anisotropic case, the off-equilibrium sim- 
ulation indicates that the spin shows the same asymptotic be- 
havior as the chirality in the isotropic case, demonstrating the 
'spin-chirality mixing' due to magnetic anisotropy. 

The authors are thankful to I. Campbell, R. Orbach, H.Takayama,| 
E.Vincent, A. P. Young, L.F.Cugliandolo, M.Ocio, H.Rieger, K. 
Nemoto and H.Yoshino for useful discussion. The numerical cal- 
culation was performed on the FACOM VPP500 at the super- 
computer center, ISSP, University of Tokyo and on the HITACHI 
SR-2201 at the supercomputer center. University of Tokyo. 
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